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Using 58 million J/ip events taken with the BESII detector at the Beijing Electron Positron 
Collider, a new decay mode J/ip — > 7/2 (1270)/2 (1270) — > 77r + -k~ tv + tt~ is observed for the first 
time. The branching ratio is determined to be Br(J/ip — > 7/2/2) = (9.5 ± 0.7 ± 1.6) x 10~ 4 , where 
the quoted errors are statistical and systematic, respectively. 

PACS numbers: 13.25.Gv, 14.40.Gx, 13.40.Hq 



I. INTRODUCTION 



Precise branching ratios are needed in order to better understand J/ip physics. Unfortunately, only 
about 50 to 60% of J / tp decay modes have been observed so far 1] . A sample of 58 million J / ip events has 
been accumulated with the upgraded Beijing Spectrometer (BESII). With this sample, the world's largest, 
it is possible to systematically study J/ip decays. 

For the decay J/tp 771-+ vr + , studies have been done by MARK III 00, DM2 and BES 
0,0. A large part of the final state is from J/ip — ► r yp°p° 0. MARK III also reported the observation of 
J/ip — *■ 777c — > 7/2/2 — * 77r + 7r~7r + 7r~ 3J. In this paper, evidence for the decay J/ip — > 7/2(1270)/2(1270) — * 
77r + 7r _ 7r + 7r _ is observed, and the branching ratio of J/ip — > 7/2(1270)/2(1270) is measured for the first time. 



II. BES DETECTOR 

BES is a conventional solenoidal magnet detector 00. A 12- layer vertex chamber (VTC) surround- 
ing the beam pipe provides trigger and trajectory information. A forty-layer main drift chamber (MDC), 
located radially outside the VTC, provides trajectory and energy loss (dE/dx) information for charged tracks 
over 85% of the total solid angle with a momentum resolution of a p /p = 0.0178-^1 + p 2 (p in GeV/c) and a 
dE/dx resolution for hadron tracks of ~ 8%. An array of 48 scintillation counters surrounding the MDC mea- 
sures the time-of-flight (TOF) of charged tracks with a resolution of ~ 200 ps for hadrons. Radially outside 
the TOF system is a 12 radiation length, lead-gas barrel shower counter (BSC). This measures the energies 
of electrons and photons over ~ 80% of the total solid angle with an energy resolution of <je/E = 21%/y/E 
(E in GeV). Outside the solenoidal coil, which provides a 0.4 Tesla magnetic field over the tracking volume, 
is an iron flux return that is instrumented with three double layers of proportional counters that identify 
muons of momentum greater than 0.5 GeV/c. 

A Geant3 based Monte Carlo, SIMBES, which simulates the detector response, including interactions 
of secondary particles in the detector material, is used in this analysis. Reasonable agreement between data 
and Monte Carlo simulation is observed in various channels tested, including e + e~ — > (7)e + e _ , e + e~ — > 
(7)/Lt/i, J/ip — > pp, J/ip — > p-K and ip(2S) — ► ir + ir~ J/ip, J/ip — > l + l~ . 



III. EVENT SELECTION 



First a J/ip — > 77r + 7r~7r + 7r~ sample is selected. Events are required to have four good charged tracks 
and one or more photon candidates. A good track, reconstructed from hits in the MDC, must be well fitted 
to a helix originating from the interaction point; have a polar angle, 9, with | cos#| < 0.8; and a transverse 
momentum greater than 60 MeV. 

The TOF and dE/dx information are used for particle identification. To identify pions, we define: 

XtofW = 



XdE/dx{i) 



(TTOF-i 



PH M Pmeasured F ^expected — i 
OPHMP-i 



Prob(i) = ProbixloFii) + xL/dxW>2). 
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Here i denotes tt, K and p. If the charged track has only TOF information or only dE/dx information, then 
Prob(i) is determined using that system only. To identify a pion, it is required that: Prob(ir) > Prob(K) 
and Prob(ir) > Prob{p). At least three tracks must be identified as pions. 

To reduce the number of spurious low energy photons produced by secondary hadronic interactions, 
photon candidates must have a minimum energy of 30 MeV and be outside a cone with a half-angle of 15° 
around any charged track. 

To get higher momentum resolution and to remove backgrounds, events are kinematically fitted to 
the J/ip — > 77r + 7r~7r + 7r~ hypothesis, looping over all photon candidates. The fit with the highest probability 
is selected, and the \ 2 01 the fit is required to be less than 20. 

For J/ip — > 77r + 7r~7r + 7r~ decay, a major source of background is from J/ip — > 7r°7r + 7r~7r + 7r~. To 
remove events containing a ir° when there are multiple photons, |m(7i72) — m(7r°)| > 60 MeV is required 

if P miss is in the plane of the two photons, 7I and 72, i.e. \P m iss ■ (?%i x ^72)! < 0.15. Here P m iss is 
the unit vector of the missing momentum of all charged tracks; f 7 i and f 7 2 are unit vectors in the 7I and 
72 directions, respectively; and 771(7172) is the invariant mass of 7I and 72. Two additional requirements, 
xt„+„-„+„- < xt^+„-„+„- and PiL < 0.0015 GeV 2 , are used to further remove J/ip — > 7r°7r + 7r~7r + 7r~ 
background. P t7 is the transverse momentum of the 7r + 7r~7r + 7r~ system with respect to the photon. Finally, 
the requirement |J7 m i S;j | = \E miss — P m iss\ < 0.07 GeV is used to reject events with multiple photons and 
charged kaons; here, E m i SS and P m iss are, respectively, the missing energy and missing momentum calculated 
using only the charged particles, which are assumed to be pions. 

There are other possible backgrounds such as J/ip — > ujir + n~, J/ip — ► jK^K^, and J/ip — > 
r yK^K ± n T . The J/ijj — > ujtt + it~ background is eliminated by the requirement \M„+„- n o — M w | > 40 MeV, 
where the tt° in ir° ir + tc~ ir + ir~ is associated to the missing momentum and energy determined using only 
the charged tracks. To remove the J/tp — > -fK®K® background, we require \M V + V - — M K o\ > 25 MeV 
for both n + n~ pairs. To remove the background from J/ip — ► jK^K^w^ , events are rejected if 
X 2 {J/ip -> jK ± tt t 'k + t:-) < x 2 (J/i/> -> 77r+7r"7r+7r") when \M n+n - - M K o \ < 25 MeV. 

Figure^a) shows the scatter plot of M^x^A versus M T27 r3 for surviving events with M 7r +„- 7r + v - > 
2.0 GeV, where irl (jr3) and 7r2 (tt4) are the ir + and ir~ with the higher (lower) momentum. There is a 
clear signal near (M p o, M p o). Since according to Monte Carlo study most events from J/ip — > 7/2/2 have 
7r + 7T~7r + 7r~ invariant mass greater than 2.6 GeV, the scatter plot of M n i n i versus M^irz with M^+^- v + v - > 
2.6 GeV is shown in Fig. Q^b). An obvious cluster near (Mf 2 , Mf 2 ) is seen besides the signal near (M„o, 
M p o). 

There is another possible M^+ n - versus M n +„~ combination, which is M^i^ versus M^^i- Monte 
Carlo studies of J/ip — > 7/2/2 and J/ip — ► jp°p° show that the j^fi and p°/5° signal almost always appears 
in the M^i^ versus M w27 r3 combination. Therefore we use the M^x-Ki versus M^-kZ combination to study 
J/ip — > 7/2/2 decay. The M^i t4 and M^i-kZ distributions for M„.+ 7r - 7r + 7r - > 2.6 GeV are shown in Fig. ^Jc) 
and Fig. ^d) respectively. Besides the p peak, a peak near M/ 2 (1.275 GeV) can be clearly seen in both of 
these plots. 

IV. BACKGROUND ANALYSIS 

Possible backgrounds are studied using simulated exclusive channels normalized according to Particle 
Data Group (PDG) branching ratios, as well as a 30 million event inclusive J/ip decay Monte Carlo sample 
generated with Lundcharm |9(. Simulated exclusive decays studied are J/ip — > r yp°p°, J/ip — > 77r + 7r - 7r + 7r~, 
J/ip — > 7r°7T + 7r~7r + 7r~, J/ip — > W7r + 7r~, J/^> — > w/2 — > ujit + it~ , and J/-0 — > 7r + 7r _ 7r + 7r~, where events 
are generated according to phase space and normalized by PDG yj branching ratios. After passing the 
simulated events through our selection criteria, no f% signal remains from these channels. 

Fig. [3 shows the M„+ n - spectra from data and from 30 million inclusive J/ip Monte Carlo decays 
normalized to 58 million events, where M 7r + 7r - 7r + 7r - is required to be greater than 2.6 GeV and the invariant 
mass of the other tt + tt~ pair satisfies \M n + n - — M/ 2 | < 0.185 GeV. This M T + T - spectrum is used to fit the 
J/ip — ► 7/2/2 signal in Section In Fig. Ufa), there is an f 2 signal from J/ip — > /2P ± 7r T decays from the 
inclusive J/ip Monte Carlo sample. After removing this channel, there is no f% signal left in the inclusive 
J/ip decay Monte Carlo events, as shown in Fig. Upa). So the f% signal in the M„+^- spectrum of inclusive 
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FIG. 1: (a) M n i„A versus M l2l3 for M w+n - n+n - > 2.0 GeV. (b) M^a versus M^s for M. 
(c) M n in4, and (d) M-ni-nz distributions with M„+„-„+„- > 2.6 GeV. 
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FIG. 2: The M r + T - spectra from data and inclusive J/ip decay Monte Carlo data, which are the sum of M^-i 7r 4 and 
M n 2Tv3- The dots with error bars are data, the shaded histogram is from inclusive J/ip decay Monte Carlo data, (a) 
Including J/ip —* fap^ir^. (b) J/ip — > f<2 Tf^ removed. 
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J/ip Monte Carlo decays is caused by J/ip — > f2P ± n T . 

Although J/ip — > /2/0 ± 7r =F decay is not listed in the PDG, we observe this decay in J/ip — > 
7r°7r + 7r~7T + 7r - , and it gives a background in the M^+ v - spectrum of J/ip — > 7/2/2- There is an / 2 
signal in the M„+ v - distribution of J/ip — > 7r°7r + 7r~7r + 7r~ when a mass requirement is made on 
the other -k^tt® pair. The number of fitted f% events in this M^+^- distribution is 2111, and the effi- 
ciency, obtained by a Monte Carlo simulation, for this J/ip — > f2P ± ^ selection is 0.67%. The number of 
J/tp — > /2P ± ti' =f — > 7r°7r + 7r _ 7r + 7r _ events is Nf 2P ± 7r T = 3.15 x 10 5 , which is used to determine the amount of 
fyp^TT^ background in the M v + V - spectrum of J/ip — > 7/2/2- 

V. RESULTS AND SYSTEMATIC ERRORS 

The Mtt+j,.- spectrum for events with M^*- > 2.6 GeV and with one tt + h~ pair satisfying 
\M„+„- — Mf 2 1 < 0.185 GeV is used to obtain the number of J/ip — > 7/2/2 events. The 7r + 7r~ mass spectrum 
is fitted with a p°, a / 2 (1270), and a 2nd order polynomial background plus /2P ± 7r T background using a 
binned maximum likelihood fit. The mass and width of the p° are fixed at the PDG values. The shape of 
the / 2 is determined by Monte Carlo simulation for J/tp -> 7/2/2 (M h = 1.275 GeV, T h = 0.185 GeV). 
The shape of the jip^it^ background is determined by Monte Carlo simulation of J/ijj — » fep ir^ and 
normalized to the estimated number of J/ip — > fap^if^ background events. The fitting result is shown in 
Fig. [3] The M„+„- spectrum is the sum of the M^i^i and the M„2ir3 spectra. The fitted number of /2 
component in the M n +„- spectrum is 1292 (both projections are used). 




M n+n -(GeV/c 2 ) 

FIG. 3: Fit of the M^+ v - spectrum to obtain the fa signal. The light shaded area is the 2nd order polynomial 
background. The dark shaded area is /2/0 ir^ background. The two peaks are p° and /2(1270), respectively. 

The branching ratio of J/ip — > 7/2/2 can be determined with the following formula: 



Br(J/i> - 7/2/2) = - 5-/ 

= (9.5 ±0.7) x 10~ 4 

where iV ob;5 =646(=1292/2) is the average number of survived J /if) — ► 7/2/2 events in M^ini projection and 
M^kZ projection, e = 3.71% is the efficiency determined from phase space Monte Carlo events with the 
same cuts used to select J/ip — > 7/2/2 data sample, Nj/^ = 57.7 x 10 6 is the total number of J/ip events 
collected by BESII an( A Br^^^^^- = 0.565 is the branching ratio of /2 — > tt + tt~ pj. 
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TABLE I: Summary of systematic errors 



Sources 


Systematic error (%) 


MDC tracking efficiency 


8 


PID efficiency 


3 


Photon selection 


2 


Kinematic fit 


4 


The number of J/tp events 


4.7 


/2 width 


5.1 




9.7 


Background uncertainty 


5.5 




2.8 


Total 


16.5 



Many sources of systematic errors are considered. Systematic errors associated with efficiencies, such 
as from the MDC tracking, particle identification, photon selection, and kinematic fitting, are determined 
by comparing ip(2S) and J/ip data with Monte Carlo simulation for very clean decay channels, such as 
ip(2S) -► TT+TT-J/ip. 

The MDC tracking efficiency has been measured using channels like J ftp — > AA and ip(2S) — ► 
7r + 7r _ J/ip, J/ip — > p + p + ■ It is found that the efficiency of the Monte Carlo simulation agrees with that 
of data within 1-2% per charged track. The total systematic error from the uncertainty of MDC tracking 
efficiency in our analysis is taken as 8%. The particle identification (PID) efficiency systematic error is 
calculated by comparing the efficiency of data with that of Monte Carlo. The largest difference between 
these two efficiencies is within 3%. According to the study of photon detection efficiency in J/ip — > pn [Tl| . 
SIMBES simulates the photon detection efficiency in the full energy range within 1 to 3%. Here, we take 
2% as the systematic error in the photon detection efficiency. 

The systematic error for the kinematic fit is caused by the differences between data and simulated 
data in the momenta and the error matrices of charged tracks and the energies and the directions of neutral 
tracks. To check the consistency between data and Monte Carlo simulation, two channels, J/tp — > p°ir° 
and J/ip — > AA, are analyzed for 2-prong and 4-prong events respectively. The systematic error for 4-prong 
events caused by kinematic fit is determined to be 4%. 

Other sources of systematic errors come from the uncertainty of Nji^, the f 2 width, the M 7T + 7T - 7T +^- 
requirement, and the uncertainty of the background. The total J/ip number is (57.7 ±2.7) x 10 6 determined 
from 4-prong events The width of the f 2 is 185.li|g MeV Q; 185.1 MeV is used in the Monte Carlo 

simulation. To estimate the systematic error caused by the width of f 2 , 188.5 MeV and 182.5 MeV are also 
tried. After refitting the 7r + 7r~ mass spectrum of Fig. [21 the relative differences of the branching ratios 
determined by using these two widths compared with the original one are 5.1% and -3.5% respectively. The 
systematic error caused by the f 2 width is taken as 5.1%. 

The requirement M n + 7r - 7I + 7I - > 2.6 GeV is used to select events for the determination of Br(J/ip — > 
7/2/2)- The Afjr+jr-Tr+jr- spectra of real data and Monte Carlo events are somewhat different. To estimate 
the systematic error caused by this requirement, we calculate Br (J/ip — > 7/2/2) with M„+ v - v + v - > 2.7 GeV 
and compare the result with that obtained with M 7r + 7r - 7r + 7r - > 2.6 GeV. The relative difference is 9.7%, 
and this is taken as one of the systematic errors. 

The background systematic error is caused by the uncertainty in the J/ip — » f 2 p tt^ background. 
This uncertainty is determined by varying the selection criteria used to obtain the number of J/ip — > f 2 p ± ^ 
events. The systematic error caused by the uncertainty of J/ip — > f 2 p ±r K T background is 5.5%. 

Table H] lists all the systematic errors, as well as the total systematic error of 16.5%, obtained by 
adding all contributions in quadrature. The branching ratio of J/ip — > 7/2(1270)/2(1270) is 

Br(J/tp -> 7 / 2 (1270)/ 2 (1270)) = (9.5 ± 0.7 ± 1.6) x 1(T 4 , 

where the first error is statistical and the second systematic. This branching ratio includes contributions 
from intermediate states such as J/ip — > jr] c . 
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VI. SUMMARY 

A new decay mode J/ip — > 7/2(1270)/2(1270) — ► 77r + 7r~7r + 7r~ is observed, and its branching ratio 
is measured to be Br(J/ip — > 7/2/2) = (9.5 ± 0.7 ± 1.6) x 10~ 4 . This result will be helpful to understand 
the complex J/tp — > -fir + ir~Tr + Tr~ decay. 
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